Non-degenerate shell-model effective interactions from 
the Okamoto- Suzuki and Krenciglowa-Kuo iteration methods 
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We present calculations of shell-model effective interactions for both degenerate and non- 
degenerate model spaces using the Krenciglowa-Kuo (KK) and the extended Krenciglowa-Kuo 
iteration method recently developed by Okamoto, Suzuki et al. (EKKO). The starting point is 
the low-momentum nucleon-nucleon interaction Vi ow -k obtained from the N 3 LO chiral two-nucleon 
interaction. The model spaces spanned by the sd and sdpf shells are both considered. With a 
solvable model, we show that both the KK and EKKO methods are convenient for deriving the 
effective interactions for non-degenerate model spaces. The EKKO method is especially desirable 
in this situation since the vertex function Z-bax. employed therein is well behaved while the cor- 
responding vertex function Q-box employed in the Lee-Suzuki (LS) and KK methods may have 
singularities. The converged shell-model effective interactions given by the EKKO and KK methods 
are equivalent, although the former method is considerably more efficient. The degenerate sd-shell 
effective interactions given by the LS method are practically identical to those from the EKKO and 
KK methods. Results of the sd one-shell and sdpf two-shell calculations for 18 0, 18 F, 19 and 
19 F using the EKKO effective interactions are compared, and the importance of the shell-model 
three-nucleon forces is discussed. 

PACS numbers: 21.60.Cs, 21.30.-x,21.10.-k 



I. INTRODUCTION 

The nuclear shell model has provided a very successful 
framework for describing the properties of a wide range 
of nuclei. This framework is basically an effective theory 
[H-Q, corresponding to reducing the full-space nuclear 
many-body problem to a model-space one with effective 
Hamiltonian PH e ffP=PH P+PV ef fP, where H is the 
single-particle (s.p.) Hamiltonian and P represents the 
projection operator for the model space which is usually 
chosen to be a small shell-model space such as the sd 
shell outside of an 16 closed core. The effective interac- 
tion V e f f plays a central role in this nuclear shell model 
approach, and its choice and/or determination have been 
extensively studied, see e.g. As discussed in these 

references, Veff may be determined using either an em- 
pirical approach where it is required to reproduce selected 
experimental data or a microscopic one where V e tf is 
derived from realistic nucleon-nucleon (NN) interactions 
using many-body methods. The folded-diagram theory 
yj-yj is a commonly used such method for the latter. 
Briefly speaking, in this theory V e f / is given as a folded- 
diagram series UHl, @ 



Q / Q-Q 



Q / Q- 



V eff = Q-Q I Q 

(1) 

where Q represents a so-called Q-box, which may be writ- 
ten as 



QH = [PVP + PVQ 



1 



u - QHQ 



QVP}! 



(2) 



Here V represents the NN interaction and ui is the so- 
called starting energy which will be explained later (sec- 
tion II) . Thus from the NN interaction V we can in prin- 
ciple calculate the Q-box and thereby the effective inter- 
action V e ff. Note that we use Q, without hat, to denote 
the Q-space projection operator. (P + Q = 1.) Note also 
that the Q-box is an irreducible vertex function where 
the intermediate states between any two vertices must 
belong to the Q space. As indicated by the subscript L 
in eq. ([2]), the Q~box contains valence linked diagrams 
only, such as the 1st- and 2nd-order Q-box diagrams for 
18 and 18 F shown in Fig. 1. The Q'-box of eq. (P) is 
defined as (Q — PVP) , namely Q begins with diagrams 
2nd-order in V. The above folded-diagram formalism has 
been employed in microscopic derivations of shell model 
effective interactions for a wide range of nuclei @, HJ ■ 

In the present work we would like to explore an ex- 
tension of the well-known methods for computing the 
folded diagram series in eq. ©. The Lee-Suzuki (LS) 0- 
9] iteration scheme has been commonly used in previous 
microscopic calculations of shell model effective interac- 
tions @,|3j]. Here we would like to employ two different 
methods, the Krenciglowa-Kuo (KK) iteration method 
[ToL [ll| and the newly developed extended Krenciglowa- 
Kuo iteration method of Okamoto, Suzuki, Kumagai and 
Fujii (EKKO)[l2j|. mainly for the purpose of calculating 
the shell-model efective interactions for non-degenerate 
model spaces. As we shall discuss later, it is not conve- 
nient to use the LS method for calculating the effective 
interactions for non-degenerate model spaces such as the 
sdpf two-shell case, while both the EKKO and the KK 



methods can be conveniently applied in this situation. 
The EKKO method has an additional advantage. When 
the P- and Q-space are not adequately separated from 
each other, the Q-box employed in the KK method may 
have singularities, causing difficulty for its iterative so- 
lution. An essential and interesting difference between 
the EKKO and KK methods is that the EKKO method 
employs the vertex function Z-hox. (to be defined in sec- 
tion II) while in the latter the vertex function Q-box 
is used. This simple replacement (of Q by Z) has an 
important advantage in circumventing the singularities 
mentioned above. As we shall discuss later, both the 
EKKO and KK methods may provide a suitable frame- 
work for calculating shell-model effective interactions for 
large non-degenerate model spaces which may be needed 
for describing exotic nuclei with large neutron excess. 

The organization of the present paper is as follows. 
In section II we shall describe a non-degenerate version 
of the EKKO metho d Jl2l| and how we apply it and 
the KK method [l^, lllj to shell-model effective inter- 
actions. A comparison of these two methods with the 
LS scheme 0, @| will be made. Our results will be pre- 
sented and discussed in section III. We shall first per- 
form a sequence of model calculations comparing the 
EKKO, KK and LS iteration methods for both degen- 
erate and non-degenerate model spaces. The use of the 
EKKO and KK methods in calculating the effective in- 
teractions for non-degenerate model spaces will be em- 
phasized. Starting from the Viow—k interaction 
fll | derived from the chiral N 3 LO potential pjj, the LS, 
KK and EKKO methods will all be used to calculate 
the degenerate sd one-shell effective interactions, a main 
purpose being to check if the results given by the com- 
monly used LS method agree with the KK and EKKO 
ones. The KK and EKKO methods will then be employed 
to calcualte the non-degenerate sdpf two-shell effective 
interactions. The matrix elements of the above degen- 
erate and non-degenerate interactions will be compared. 
The low-energy spectra of 18 0, 18 F, 19 and 19 F given 
by these interactions will be discussed. A summary and 
conclusion will be presented in section IV. 



II. FORMALISM 

In this section, we shall describe and discuss the KK 
ToL 11 1 and the EKKO 12 1 iteration methods and their 
application to microscopic calculations of shell-model ef- 
fective interactions. These methods, to our knowledge, 
have not yet been employed in such calculations. Let us 
begin with a brief review of the LS and KK iteration 
methods. Consider first the degenerate LS method where 
the model space is degenerate, namely PH P — Wo, Wq 
being a constant. In terms of the Q-box of Eq.(2), the ef- 
fective interactions R n are calculated iteratively by 0, H| 



'WWn 
\ZVW-' 






FIG. 1: Low-order diagrams constituting the Q-box. 



Ri = 



1-Q 



-Q(W ), 



and for n > 2 

Rn = 



1 Ql Em=2 Q m Ylk=n-m+l 



(4) 



Q(W ), (5) 



where Q rn is proportional to the mth derivative of Q: 



1 d m Q 



il du r 



(6) 



The effective interaction is given by the converged R n , 
namely 



^e// — Rn+1 = Rn- 



(7) 



There is a practical difficulty for the above iteration 
method. In actual shell-model calculations, it is usually 
not possible to calculate the vertex function Q-box ex- 
actly; thus it is a common practice to evaluate it with 
some low-order approximation and calculate the deriva- 
tives Q m numerically. At higher orders this becomes in- 
creasingly difficult, and therefore such calculations are 
usually limited to low orders in the iteration. However, 
as we shall demonstrate in section Ilia, low-order LS iter- 
ations are often not accurate when the P- and Q-spaces 
are strongly coupled. 

The above degenerate LS iteration method can be gen- 
eralized to a non-degenrate one [jj, namely PHqP being 
non-degenerate. In this situation we need not only the 
Q-box of eq. @ but also a generalized Q-box defined by 



Qn{eie 2 ■ •■£„+!) = (-l) n [PVQgfgf ■ --g^QVPU 
with 



e t - QHQ 



(8) 
(9) 



Ri = Q(W ), 



(3) 



This generalized Q-box is defined for n > 1, and ti is de- 
fined by PHoPfa = ei<pi where P = 52 m<d \<p m ){<f>m\- 



3 



The dimension of the P-space is labelled d. Note 
that only valence-linked diagrams are retained in 
On( e i • • • as indicated by the subscript L. 

With the above definitions, the effective interaction 
given by the non-degenerate LS iteration method is given 
as [| 



Ri = ^Q(e Q )P Q , 



R 2 = £ 



1 - ]T Qi (e afJ ) Pp 



* = E 

a 

xQ(e Q )P ( 



Pi 



[PH Q P 



V { %]r {i) = E [l 



x 



(E m (uj) - H )xm = Q(u)Xm, w = E m (u). 
To start the iteration, we use 



(i) 



Q(wo) 



where wo is a starting energy chosen to be close to PHqP. 
The converged KK effective interaction is given by V e f / = 

y^t^ ~ K//- When convergent, the resultant Veff is 
independent of ujq, as it is the states with maximum P- 
space overlaps which are selected by the KK method [l(| ■ 
We shall discuss this feature later in section III using a 
solvable model. The above non-degenerate KK method 
is numerically more convenient than the non-degenerate 
LS method. 

The diagrams of Fig. 1 have both one-body (dl,d2,d3) 
and two-body (d4,d5,d6,d7) diagrams. When we calcu- 
late nuclei with two valence nucleons such as 18 and 
i 18 F, all seven diagrams are included in the Q-box. But 
for nuclei with one valence nucleon such as 17 O, we deal 
with the 1-body S'-box which is approximated by the sum 
of diagrams dl, d2 and d3. The 1-body effective interac- 
tion is given by a similar KK iteration 



(10) 



with e a p = (e Q + ep), e a p y = (e a + ep + e 7 ), • • •. In 
the above equations P m = \4>m){4>m\- When convergent, 
we have V e ff = Rn+i = Rn- The above non-degenerate 
LS method can be used to calculate the effective inter- 
actions for e.g., the non-degenerate Op-shell jl3[ or the 
two-shell sdpf model space. But this method is rather 
complicated for computations, and this has hindered its 
application to microscopic calculations of shell-model ef- 
fective interactions. 

We now describe some details of the non-degenerate 
KK and EKKO methods. The KK iteration method was 
originally developed for model spaces which are degen- 
erate [lOj. A non-degenerate KK iteration method was 
later formulated ll| , with the effective interaction given 
by the following iteration methods. Let the effective in- 
teraction for the ith iteration be V^ff and the correspond- 
ing eigenvalues E and eigenfunctions x be- given by 



(11) 



Here Xm is the P-space projection of the full-space eigen- 
function namely Xm = P^m- The effective interac- 
tion for the next iteration is then 

V e}T ] =^ PH o P + Q( E ^\^)(xt ) \~PH P, (12) 

m 

where the bi-orthogonal states are defined by 



(13) 



Note that in the above PHqP is non-degenerate. The 
converged eigenvalue E m and eigenfunction x-m satisfy 
the P-space self-consistent condition 



(14) 



(15) 



S eff =J2i PH ^ P + § ( E ^\^)^\- PH R ( 16 ) 
m 

Denoting its converged value as S e ft, the model-space 
s.p. energy e^ f is given by P m (H + S ef f)P m . By adding 
and then subtracting S e ff, we can rewrite Eq.(14) as 



Hn ff )Xr 



[Q{oj) - S eff }x m , w = E m (w), 

(17) 

with Hf) = Hq + S e ff. In most shell model calcula- 
tions [Jjrj, one often uses the experimental s.p. energies. 
This treatment for the s.p. energies is in line with the 
above subtraction procedure, as P(Hq + S e ff)P repre- 
sents the physical s.p. energy which in principle can be 
extracted from experiments. In the present work we shall 
use the experimental s.p. energies for the model-space or- 
bits together with the V e }f derived from (Q — S e ff). A 
similar subtraction procedure has also been employed in 
the LS calculations [2j, |3( where the all-order sum of the 
one-body diagrams was subtracted from the calculation 
of the effective interaction. 

In several aspects, the above KK method provides a 
more desirable framework for effective interaction cal- 
culations than the commonly used LS method. The 
KK method is more convenient for non-degenerate model 
spaces than the LS method, and the KK method does not 
require the calculation of high-order derivatives of the 
Q-box, which may be necessary in a converged LS cal- 
culation. The KK method has, however, a shortcoming 
when applied to calculations with extended model space 
such as the two-shell sdpf one. For example, certain 
2nd-order diagrams for this case may diverge, resulting 
in an ill-defined Q-box. It is remarkable that these po- 
tential divergences can be circumvented by the recently 
proposed EKKO method of Okamoto et al. [12]. In this 
method, the vertex function Z-box is employed. It is 
related to the Q-box by 



1 

i-QiH 



[<3(w)-&(w)P(w-.Ho)P], (18) 
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where Qi is the first-order derivative of the Q-box. The 
Z-bax. considered by Okamoto et al. [12] is for degener- 
ate model spaces (PH P — W ), while we consider here 
a more general case with non-degenerate PH P. An im- 
portant property of the above Z-box is that it is finite 
when the Q-box is singular (has poles). Note that Z(ui) 
satisfies 



Z(w)Xm = Q(u)Xm at uj = E m (uj). 



(19) 



The iteration method for determining the effective inter- 
action from the Z-box is quite similar to that for the 
Q-box. Suppose the effective interaction for the ith iter- 
ation is V^_ z . The corresponding eigenfunction \ and 
eigenvalues E z are determined by 

[PH P + V$ f _ z ] X ® = E^xtt- (20) 
The effective interaction for the next iteration is 

V Sf-z = }2l PH o P + Z{E^)]\x^){m ~ PHoP, 
m 

(21) 

Although Q{uS) and Z(ui) are generally different, it is in- 
teresting that the converged eigenvalues E m of P(H + 
V e ff)P and the corresponding ones Ef n of P(H + 
V e ff_z)P are both exact eigenvalues of the full-space 
Hamiltonian H — Hq + V, which can be seen from eqs. 
(14), (18) and (19). Note, however, the KK and EKKO 
methods may reproduce different eigenvalues of the full- 
space Hamiltonian H. This aspect together with some 
other comparisons of these methods will be discussed in 
section Ilia, using a simple solvable model. 

For the degenerate case, Okamoto et al. [IH have 



shown that 



=0 at the self-consistent point u> 



E^(u>). As outlined below, we have found that this result 
also holds for the case of non-degenerate PHqP. From 
eq. (18), we have 



dZ 
du 



1 



1-Qi l-Qi 
1 — (Ji 



Then from eqs. (14) and (19) we have 



Q-Qx(u-H )] (22) 



^dZ{uj) 
dio 



]ul=EZ \Xm) = 0, 



and 



, dEg(u) 
1 du 1 



= 0. 



(23) 



(24) 



This is a useful result; it states that at any self-consistent 
point the eigenvalues E^(uS) varies 'flatly' with w, a fea- 
ture certainly helpful to iterative calculations. In section 
III, we shall check this feature numerically. 



III. RESULTS AND DISCUSSION 



Ilia. Model calculations comparing the LS, KK and 
EKKO methods 



In this section we shall study the above iteration meth- 
ods by way of a simple matrix model, similar to the one 
employed in [l2| ■ We consider a 4-dimensional matrix 
Hamiltonian H = Hq + H% where 



and 



PH Q P 



Ho 



£pl 
£p2 



PH a P 
QH Q 



; QH Q Q 



£ql 
£q2 



(25) 



(26) 



The interaction Hamiltonian has a strength parameter x, 
namely 



Hi 



PHxP PPlxQ 
QHiP QHiQ 



(27) 



with 



bx 
5a; lOx 



PHxP = 

PH X Q = ()UyP = 

QHiQ = 



8a; 
8x 



-5x 

x 



x 
5x 



(28) 



As discussed in section II, both the KK and EKKO it- 
eration methods are rather convenient for non-degenerate 
model spaces. We would like to check this feature by 
carrying out some calculations using the above model. 
We consider two unperturbed Hamiltonians, given by 
(e P i, £ P 2, £ g i, e g2 )=(0,6,4,9) and (0,0,4,9). The PH P 
parts of them are, respectively, non-degenerate and de- 
generate. 

In the first three entries of Table I, some results for 
the PH o P—(0fi) case are presented. Here E n are the ex- 
act eigenvalues of the full Hamiltonian, with their model- 
space overlaps denoted by (n\P\n). The entries Ekk and 
Eekko are the eigenvalues generated respectively by the 
KK and EKKO iteration methods. Not only is the above 
PHqP non-degenerate but its spectrum intersects that 
of QHoQ. One would expect that this PHqP may cause 
difficulty for the above iteration methods. But as indi- 
cated in Table I, both the non-degenerate KK and the 
non-degenerate EKKO iteration methods work remark- 
ably well. Note that the interaction used here is rather 
strong (x=0.6), and both methods still work well, con- 
verging to values of E n which are quite far from PHqP. 

Some results for the degenerate case of PHqP=(0,Q) 
are listed in the last two entries of Table I. Here we 
have performed calculations using the degenerate LS 
method through 5th order iteration (i.e. in Eq.(7) we use 
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V e ff = -R5). As shown, the results so obtained are not 
in good agreement with the exact results. This suggests 
that low-order LS iteration method may often be inad- 
equate, and one needs higher-order iterations to obtain 
accurate results. 

It is known that the KK iteration method converges to 
the states with maximum P-space overlaps floT ], while the 
LS method converges to the states of lowest energies 0, 
8] . We have found that for many cases the EKKO method 
also converges to states of maximum P-space overlaps. 
As listed in the third entry of Table I, both the KK and 
EKKO methods converge to states of energies E n = -3.51 
and 14.53 whose P-space probabilities are relatively 0.70 
and 0.86. We have also found that the EKKO and KK 
iteration methods can converge to different states. An 
example is the result shown in the last part of Table I, 
where the EKKO method converges to states of energy 
(P-space probability) -1.45 (0.87) and 0.91 (0.46), while 
the states of maximum P-space overlap are those with 
energies -1.45 and 5.25. Note that for this case the EKKO 
method is clearly more accurate than the KK method. 

We have noticed that for a number of cases the EKKO 
method converges well but not so for the KK method. 
This is largely because these two methods treat the sin- 
gularities of the Q-box differently. To see this, let us 
perform a graphical solution for the x=0.60 case of Table 
I. Using the parameters of this case, we calculate and plot 
in Fig. 3 both E^(cj) and E z z (uj) which are respectively 
the eigenvalues of P[H +Q(u)]P and P[H +Z(u)]P. As 
discussed in section II, they have identical self-consistent 
solutions, namely 10 = E®(uj) = E z n {uj) = E m where E m 
is the eigenvalue of the full-space Hamiltonian. Recall 
that Q and Z are given respectively by eqs. (2) and (18). 
As shown in the figure, the curves of E® and E z do have 
identical self-consistent solutions as marked by the com- 
mon intersection points E\, P2, £3 and E±. Note that 
the above two curves are distinctly different from each 
other, particularly in the vicinity of the poles (marked 
by the vertical lines through Pi and P2) of the Q-box. 
There E^(u) is discontinuous, diverging oppositely be- 
fore and after the pole, while E z (lu) remains continuous 
throughout. This clearly helps the convergence of the 
Z-box iteration method: The Z-box iteration proceeds 
along a continuous E z (w) curve, while the Q-box itera- 
tion often does not converge as it may bounce back and 
forth across the discontinuity. 

As seen from eq. (18), the Z-box method has 'false' 
solutions at E z (w) = oj = [i q = F q , where \i q are poles 
of the Q-box. These solutions are marked in Fig. 3 as 
Pi and P2. These false solutions can be readily recog- 
nized and discarded. As given in Eq.(24), we have at 

self-consistent points = 0. As shown in Fig. 3, the 
slopes of E z do satisfy the above condition at the self- 
consistent points E\ to P4, but not so at the false points 
Pi and P 2 . 



TABLE I: Results of model calculations using the LS, KK and 
EKKO iteration methods. See text for other explanations. 



H = (0,6,4,9) 


x=0.10 








E,, 


-0.110705 


3.328164 


7.203243 


8.579299 


(n\P\n) 


0.991108 


0.044558 


0.953368 


0.010966 


Ekk 


-0.110705 


7.203243 






Eekko 


-0.110705 


7.203242 








x =0.30 








E n 


-0.974184 


1.808716 


8.080888 


10.084581 


(n\P\n) 


0.906155 


0.111257 


0.113523 


0.869066 


Ekk 


-0.974185 


10.084580 






Eekko 


-0.974185 


10.084581 








x =0.60 








E n 


-3.510377 


-0.286183 


8.265089 


14.531471 


(n\P\n) 


0.709531 


0.189318 


0.232422 


0.868729 


Ekk 


-3.510377 


14.531472 






Eekko 


-3.510363 


14.531472 






F =(0,0,4,9) 


x= 0.10 








-En 


-0.296201 


0.982861 


3.736229 


8.577111 


(n\P\n) 


0.985416 


0.922993 


0.082811 


0.008780 


Els 


-0.314591 


0.872176 






Ekk 


-0.296201 


0.982860 






Eekko 


-0.296201 


0.982860 








x= 0.30 








E n 


-1.448782 


0.906504 


5.254129 


8.288149 


(n\P\n) 


0.868986 


0.460296 


0.566014 


0.104704 


Els 


-1.620734 


0.339561 






Ekk 


-1.125579 


5.515468 






Eekko 


-1.448782 


0.906504 
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FIG. 2: Graphical solutions for the Q- and Z-box self- 
consistent equations. See text for other explanations. 
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FIG. 3: Energy spectrum of 1 O obtained from the sd-shell 
calculations using the LS, KK and EKKO methods. See text 
for more explanations. 



Illb. The sd and sdpf shell model effective 
interactions 



In this subsection, we shall calculate the effective in- 
teractions for both the degenerate sd one-shell and the 
non-degenerate sdpf two-shell cases. Before presenting 
our results, let us first describe some details of our calcu- 
lations. The LS, KK and EKKO methods as described 
in section II will be employed. We first compute the low- 
momentum nucleon-nucleon interaction V^_ k [1, 14- 16 
starting from the chiral N 3 LO two-body potential 17 
at a decimation scale of A = 2.1 frn -1 . At this cutoff 
scale, the low-momentum interactions derived from dif- 
ferent NN potentials (r7l - [20| are remarkably close to each 
other, leading to a nearly unique low-momentum inter- 
acton 16]. The effect of the leading-order chiral three- 
nucleon force on shell model effective interactions has 
been studied, and our results will be reported in a sep- 
arate publication [2l|. The above V^_ k interaction is 

then used in calculating the Q-box diagrams as shown 
in Fig. 1. In calculating these diagrams, the hole or- 
bits are summed over the OsOp shells and particle orbits 
over the OdlslpOf shells. The active spaces (P-space) 
used for the one-shell and two-shell calculations are re- 
spectively the three orbits in the sd shell and the seven 
orbits in the sdpf shells. The experimental s.p. energies 
of (0.0, 5.08, 0.87) MeV have been used, respectively, for 
the (0^5/2; 0(^3/2, lsi/2) orbits [22|- We have employed 
the shell-model s.p. wave functions and energies with the 
harmonic oscillator constant of huj— 14 MeV. 

In the following, let us first report our results for the 
degenerate sd one-shell calculations. Here by degenerate 
we mean that the unperturbed s.p. energy levels for the 
((W5/2, 0c?3/2, lsi/2) orbits are degenerate. Our purpose 
here is mainly to compare the results given by the KK 
and EKKO method with those given by the commonly 
used degenerate LS method 0, H| ■ Our results are pre- 



FIG. 4: Same as Fig. 3 except for 18 F. 



scnted in Figs. 3 and 4. As discussed in section II and 
illustrated in Table I, the EKKO and KK methods may 
converge to different states. In the present calculations 
they actually converge to the same states, as seen in Figs. 
3 and 4. Our LS calculations are carried out using a low- 
order approximation, namely we take V e f / = R5 (see eq. 
(7)). As shown in the figures, the LS, KK and EKKO 
results for both 18 and 18 F are in fact nearly identical 
to each other. A comparison of our results with exper- 
iments Q is also presented in the figures. The agree- 
ment between our calculated energy levels with experi- 
ments is moderately satisfactory for 18 0, but for 18 F the 
calculated lowest (l + ,3 + ,5 + ) states, though of correct 
ordering, are all significantly higher than the experimen- 
tal values. As discussed in sections II and Ilia, the LS 
method is known to converge to the states of the lowest 
energies, while the KK method to the states of maxi- 
mum P-space overlaps. Thus the good agreement shown 
in Figs. 3 and 4 is an important indication that the states 
reproduced by the model-space LS, EKKO and KK ef- 
fective interaction are likely those of the lowest energies 
as well as maximum model-space overlaps. In Fig. 5, 
we compare the entire sd-shell matrix elements given by 
the LS and EKKO methods; it is remarkable that every 
individual LS matrix element is practically equal to the 
corresponding EKKO one. (The matrix elements given 
by EKKO and KK are also nearly identical.) Recall that 
our LS matrix elements were obtained with a low-order 
i?5 iteration, which indicates the rapid convergence of the 
iteration scheme in the case of srf-shcll effective interac- 
tions. 

As discussed in sections II and Ilia, the EKKO or KK 
methods are convenient for deriving the effective inter- 
actions for non-degenerate model spaces. In the follow- 
ing let us first apply these methods to a relatively sim- 
ple case, namely the non-degenerate sd one-shell effec- 
tive interactions. The shell-model s.p. energies for the 
sd shell are degenerate, but the experimental ones are 
not. It may be of interest to employ the experimental 
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FIG. 5: Comparison of the sd-shell matrix elements given by 
the LS (LS-eff) and EKKO (KK-eff) methods. The diagonal 
line indicates the equivalence of the LS and EKKO matrix 
elements. See text for more explanations. 



FIG. 7: Comparison of the matrix elements of the sd one-shell 
(1-eff ) and sdpf two-shell (2-eff) EKKO effective interactions. 
Only the matrix elements within the sd shell are shown. See 
text for more explanations. 
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FIG. 6: Comparison of the low-lying states of 19 F calculated 
with the degenerate (DEG) and non-degenerate (NON-DEG) 
sd one-shell EKKO effective interactions. 



s.p. energies as the unperturbed s.p. energies for the sd 
shell Thus here we employ the same unperturbed s.p. 
spectrum as the previous degenerate sd case except that 
the 0g?3/2 and Is 1/2 orbits are shifted upward by, respec- 
tively, 5.08 and 0.87 MeV relative to the Ocis/2, mimicking 
the experimental s.p. energies. We have found that the 
non-degenerate and degenerate sd EKKO effective inter- 
actions do not differ significantly from each other. To 
illustrate, we compare in Fig. 6 the spectra of 19 F for 
both the degenerate and non-degenerate sd-shell calcu- 
lations; they agree with each other rather well. Slightly 
better agreements between the spectra of 18 0, 18 F and 
19 O calculated with these two interactions are also ob- 
served. In short, the one-shell sd effective interactions 
given by the degenerate and non-degenerate choices for 
the unperturbed sd s.p. energies are nearly the same, and 
it is adequate in this case to just use the former choice. 
We now report our results for the non-degenerate sdpf 



FIG. 8: Low-lying states of F calculated with the sd one- 
shell (lshell) and sdpf two-shell (2shell) EKKO interactions. 
Two choices for the vertex function Q-box are employed. For 
the left-most two columns only the first-order Q-box diagrams 
are included, while for the middle two columns we include 
both the first- and second-order diagrams. 



two-shell effective interactions. So far our calculations 
have all been carried out using the sd one-shell model 
space. For certain nuclei such as those with a large neu- 
tron excess, a larger model space such as the sdpf one 
may be needed. It will be convenient to describe our 
sdpf calculations by way of an example, namely 18 0. 
Consider the (T = 1, J = 0) states of this nucleus. In the 
sd one-shell case, the model space is spanned by three 
basis states \j 2 ,T = 1, J = 0), where j = d 5 /2,d 3 / 2 , or 
si/2- For the sdpf case, the model space is enlarged, 
having four additional single-particle states labeled by 
J = (P3/2)Pi/2>/7/2)/5/2)- The whole model space is 
now non-degenerate as the pf orbits are one shell above 
the sd ones. Similar to the sd one-shell case, we have 
calculated the sdpf effective interactions using both the 
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EKKO and KK methods. For the sd shell-model calcu- 
lations, we have employed the experimental s.p. energies 
for the three sd-shell orbits as mentioned earlier. 

For the sdpf shell-model calculations, we need in addi- 
tion the experimental s.p. energies for the four pf orbits. 
Their values are, however, not well known. In the present 
calculation we have placed them all at a separation of one 
fibj (14 MeV) above the d§ n level. (As to be reported 
later (Fig. 11), we shall also use a smaller value for the 
above separation.) In our calculations we have consid- 
ered two choices for the unperturbed s.p. energies of the 
sd model space, a degenerate one and a non-degenerate 
one (in which the 0c? 3 / 2 and lsi/2 orbits are shifted higher 
in energy compared to the 0d 5 / 2 orbit as described ear- 
lier). We have found that the results are rather similar, 
and in the following discussion we report only the calcu- 
lations for the degenerate sd-shell choice. 

In Fig. 7 we compare the matrix elements of the sdpf 
two-shell EKKO interaction with those of the sd one- 
shell case. Only the matrix elements within the sd shell 
are shown. We find that the magnitudes of the two-shell 
matrix elements are generally weaker than the one-shell 
matrix elements, some differing by as much as 1 MeV. 
Despite these differences, the resulting spectra for 18 
and 18 F given by the one-shell and two-shell calculations 
are nearly equivalent to each other, as illustrated by the 
two middle columns of Fig. 8. In our present and sub- 
sequent calculations we employ a low-order Q-box con- 
sisting of the 1st- and 2nd-order diagrams of Fig. 1. It is 
instructive, however, to compare the one- and two-shell 
calculations when only the leading-order approximation 
to the Q-box is retained. In the two leftmost columns 
of Fig. 8 we show the spectra of 18 F when only the 1st- 
order Q-box diagrams are included. We note that in this 
case (in contrast to the second-order Q-box calculation) 
the resulting sd '1 Shell' and sdpf '2 Shell' spectra are 
largely different. 

This important observation can be explained as fol- 
lows. In general, the different sd and sdpf model spaces 
will result in different effective interactions and associ- 
ated effective Hamiltonians, which we denote by -H^// 
and -fffyj- But since the sd model space is a subspace 
of sdpf model space, and H^, » should in princi- 

ple have common eigenvalues. We would like to check 
that this requirement is satisfied by our present calcula- 
tions. As indicated by the two middle columns of Fig. 
8, we see that indeed this is the case. In fact among the 
many states given by the sdpf calculations, it is the ones 
with the maximum sd-space overlaps which agree with 
the results given by the s<i-shell calculation, which is a 
physically desirable result. Because of the difference in 
the model spaces, the renormalization effects for the sd 
and sdpf effective interactions are different. To bridge 
these differences, we need to include at least the 2nd- 
order Q-box diagrams (note that the allowed intermedi- 
ate states of the 2nd-order Q-box diagrams are model- 
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FIG. 9: Low-lying states of 19 calculated with the sd one- 
shell (1 Shell) and sdpf two-shell (2 Shell) EKKO interac- 
tions. The experimental results are from [221 ]. 
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FIG. 10: Same as Fig. 10 except for 19 F. 



space dependent). These diagrams are not included in 
the above lst-order calculations, and consequently the 
sd and sdpf results are different as shown by the two 
leftmost columns of Fig. 8, which is a strong evidence for 
the importance of the model-space-dependent renormal- 
ization effects. The construction of model-space effective 
interactions is in many ways similar to the construction 
of an effective theory from a renormalization group evolu- 
tion. In such cases, one would expect that despite the dif- 
ferent effective Hamiltonians, the same low-energy phys- 
ical observables would be reproduced. The shell model 
effective interactions in the present work have not been 
computed exactly (that is, including high-order diagrams 
in the Q-box), yet it is interesting that we nevertheless 
find excellent agreement between the one- and two-shell 
calculations including Q-box diagrams only up to second 
order. 

To further study the sd one-shell and sdpf two-shell ef- 
fective interactions, we have applied them to shell model 
calculations of 19 and 19 F. The 2nd-ordcr Q-box of Fig. 



9 



> 

CD 



CD 
LU 



-10 




FIG. 11: Low-lying states of 19 and 19 F calculated using the 
sdpf two-shell EKKO interactions, with the pf experimental 
s.p. energies placed at 10 ('pf-P) and 14 MeV ('pf-II') above 
the 0d 5 /2 orbit. 
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FIG. 12: A shell-model three-nucleon force diagram. Its ex- 
ternal lines all belong to the sd shell, while its intermediate 
states between any two vertices must have at least one line 
belonging to the pf shell. 



1 is employed. Our results are displayed respectively in 
Figs. 9 and 10. It is of interest that the one-shell and 
two-shell results for 19 are nearly identical, and for 19 F 
they are also remarkably close to each other. Although 
the orderings of our calculated spectra are in fair agree- 
ments with experiments, there are significant differences 
between them. As discussed earlier, in our sdpf two- 
shell calcualtions the experimental s.p. energies for the 
four pf orbits are needed, but their values are not well 
known. So far we have chosen to place them at a separa- 
tion of 14 MeV above the 0d 5 / 2 orbit. We have repeated 
our calculations using instead a smaller separation of 10 
MeV, to investigate if the use a smaller separation may 
improve the agreements. As illustrated in Fig. 11, the 
low-lying states of lg O are hardly changed, those given 
by the separations of 14 and 10 MeV being nearly iden- 
tical. The differences between the two sets of states for 
19 F are also generally small. That these low-lying states 
are insensitive to the above separations indicates that our 
sdpf effective interactions (obtained with the inclusion of 
the 1st- and 2nd-order Q-box diagrams of Fig. 1) have a 
rather weak coupling between the sd and fp shells. 



Recall that we have employed the folded-diagram ex- 
pansion of eq. (1) to calculate the effective interaction 
V e ff. For nuclei with three valence nucleons, this expan- 
sion has both 2-body and 3-body diagrams. (By 3-body 
diagrams we mean those valence-linked diagrams with 
three incoming and three outgoing valence lines [l[.) As 
an example, the shell-model three-nucleon force diagram 
of Fig. 12 should be included in the sd one-shell effec- 
tive interaction for 19 O and 19 F. But this diagram is not 
included in our present sd one-shell calculation as we em- 
ploy only the two-body Q - box of Fig. 1. This diagram is, 
however, included in our sdpf shell model calculations 
for 19 and 19 F. Thus the good agreement between the 
one-shell and two-shell results shown in Fig. 9 is an in- 
dication that this type of shell-model three-nucleon force 
is likely to have little importance for 19 O spectra, while 
it is moderately important for 19 F as suggested by the 
small difference shown in Fig. 10. Further studies of this 
type of three-nucleon force will be useful and of interest, 
and we plan to do so in a future study. 



IV. SUMMARY AND CONCLUSION 

We have applied the iteration method of Krenciglowa 
and Kuo (KK) and that recently developed by Okamoto 
et al. (EKKO) to the microscopic derivation of the sd 
and sdpf shell-model effective interactions using the low- 
momentum nucleon-nucleon interaction derived from the 
chiral N 3 LO two-body potential. We first considered a 
solvable model and found that both methods are suit- 
able and efficient for deriving the effective interactions 
for non-degenerate model spaces, where the Lee-Suzuki 
iteration method is considerably less convenient. Even 
in the situation where the P- and Q-space unperturbed 
Hamiltonians have spectrum overlaps did the KK and 
EKKO methods perform remarkably well. The EKKO 
method has the special advantage that its vertex func- 
tion Z-box is, by construction, a continuous function of 
the energy, while the Q-t>ax function used in the LS and 
KK methods may have singularities. This feature was 
found to be particularly useful for the convergence of the 
EKKO iteration method. 

Using the Viow—k low-momentum potentialss obtained 
from above two-body interaction, we first calculated the 
degenerate sd one-shell effective interactions using the 
LS, KK and EKKO methods. The results given by KK 
and EKKO were found to be identical. It is notewor- 
thy that the LS results, calculated with a low-order (5th 
order) iteration, were also in very good agreement with 
both the KK and EKKO results, supporting the accuracy 
of the low-order LS method for calculating the degenerate 
shell-model effective interactions. 

We have calculated the non-degenerate sdpf two-shell 
effective interactions using both the EKKO and KK 
methods. Both methods gave identical results and were 
found to be suitable for such non-degenerate calcula- 
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tions, with the former being more efficient (faster con- 
verging). We have applied these interactions to compute 
the low-lying energy spectra for several nuclei with two 
and three valence nucleons above the 16 core. Since the 
sd model space is a subspace of the sdpf model space, 
we expect the effective Hamiltonians for these two spaces 
to have common eigenvalues. Indeed this was largely 
confirmed in our calculations of 18 0, 18 F, 19 and 19 F 
spectra, where it was found that the states in the sdpf 
calculations with the maximum sd-space overlap agreed 
with the results given by the sd calculations. The above 
agreement was found to be excellent for 19 O, though not 
as good for 19 F, which indicates that the shell-model 
three-nucleon force is more important in 19 F (where the 
proton- neutron interaction is involved) than in 19 O. Fur- 
ther study of this three-nucleon force should be useful 
and of much interest. 

The calculated ground state energies for the above four 
nuclei are all higher (less bound) than the corresponding 
experimental values. We are studying if the inclusion of 
the chiral three-nucleon force may give additional bind- 
ing energy [21] . In the present work we have employed 



the Q-box irreducible vertex function consisting of the 
first- and second-order diagrams only. The inclusion of 
the Kirson-Babu-Brown (KBB) all-order core polariza- 
tion diagrams in the Q-box may provide additional bind- 
ing energy [23| . We plan to carry out futher calculations 
with the inclusion of such all-order KBB diagrams. 
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under contracts DE-FG02-88ER40388 and the DFG 
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Origin and Structure of the Universe are gratefully ac- 
knowledged. 



[2 
[3; 

[4: 
[5; 
[e; 

[7. 
[8 

[9 
[10 

[11 
[12 



[13 
[14 

[15; 

[ie; 
[17; 
[is; 

[19 

[20 

[21 

[22 
[23 



M. Hjorth- Jensen, T. T. S. Kuo, and E. Osnes, Phys. 
Rep. 261, 126 (1995), and references therein. 
L. Coraggio, A. Covello, A. Gargano, N. Itaco and T.T.S. 
Kuo, Prog. Part. Nucl. Phys., 62 (2009) 135, and refer- 
ences quoted therein. 

B.A. Brown and B.H. Wildenthal, Ann. Rev. Nucl. Part. 
Sci. 38, 29(1988). 

B.A. Brown, W.A. Richter, Phys. Rev. C74, 0343150 
(2006). 

T.T.S. Kuo, S.Y. Lee and K.F. Ratcliff, Nucl. Phys. 
A176 (1971) 172. 

S.Y. Lee and K. Suzuki, Phys. lett. 91B, 173 (1980). 
K. Suzuki and S.Y. Lee, Prog. Theor. Phys. 64, 2091 
(1980). 

K. Suzuki, R. Okamoto, P.J. Ellis and T.T.S. Kuo, Nucl. 
Phys. A567 (1994) 576. 

E.M. Krenciglowa and T.T.S. Kuo, Nucl. Phys. A235, 
171 (1974). 

T.T.S. Kuo, F. Krmpotic, K. Suzuki and R. Okamoto, 
Nucl. Phys. A582, 205(1995). 

R. Okamoto, K. Suzuki, H. Kumagai and S. Fujii, 
to be published in 'Proceedings of the 10th Interna- 
tional Spring Seminar on Nuclear Physics (May 21-25, 
Sur Mare Vietri, Italy, ed. by A. Covello)'; [nucl-th] 
larXiv:1011.1994V l. 

L. Coraggio and N. Itaco, Phys. Lett. B616, 43 (2005). 
S. K. Bogner, T. T. S. Kuo and L. Coraggio, Nucl. Phys. 
A684, (2001) 432. 

S.K. Bogner, T.T.S. Kuo, L. Coraggio A. Covello, and 
N. Itaco, Phys. Rev. C 65, 051301(R) (2002). 
S. K. Bogner, T. T. S. Kuo, and A. Schwenk, Phys. Rep. 
386, 1 (2003). 

D. R. Entem, R. Machleidt, and H. Witala, Phys. Rev. 
C 65, 064005 (2002). 

R. Machleidt, Phys. Rev. C 63, 024001 (2001). 

R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. 

Rev. C 51, 38 (1995). 

V. G. J. Stoks, R. A. M. Klomp, C. P. F. Terheggen, and 

J. J. de Swart, Phys. Rev. C 49, 2950 (1994). 

H. Dong, T.T.S. Kuo and J. W. Holt, preprint (May 2011, 

to be submitted to arXiv and PRC). 

http://www.nndc.bnl.gov/chart/ 

J.D. Holt, J. W. Holt, T. T. S. Kuo, G. E. Brown and S. 
K. Bogner, Phys. Rev. C72, 041304(R) (2005). 



[1] T. T. S. Kuo and E. Osnes, Lecture Notes in Physics 
( Springer- Verlag, New York, 1990), Vol. 364. 



